The role of the band gaps in reconstruction of polar surfaces and interfaces 
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Density functional theory applied to a simple ionic material, MgO, is used as a model system to 
clarify several aspects of electronic driven mechanism to compensate for the diverging electrostatic 
potential in the polar structures. We demonstrate that in the stoichiometric free standing films, the 
electronic reconstruction is limited by the band gap. This produces a residual electric field in the 
bulk of the sample which is extremely sensitive to tiny deviations in electron transfer between two 
surfaces of the slab. In heterostructures, the band gap is replaced by a new effective energy scale 
set by the band alignment of its component. This changes the reconstruction pathways so that the 
electronic mechanism can benefit from the smallest energy scale possible. 



Beyond a doubt, the electrostatic potential plays a fun- 
damental role in the surface and bulk properties of vir- 
tually any ionic material. It has been long recognized as 
the driving force behind the so-called polar catastrophe 
that controls the physics and chemistry of polar surfaces 
and interfaces'^. A prominent example is a discovery 
of a 2D electron gas at the interface between LaAlOs 
and SrTiOs^. In the bulk, each one of these perovskites 
is an insulator with large conductivity gap of 5.6 and 
3.25eV, respectively. However, when the two are com- 
bined in a heterostructure such that the (001) plane be- 
comes an interface, a new metallic state is created^. Sev- 
eral proposals have been put forward to explain the ef- 
fect. It is argued that the large sheet charge density 
could be a result of Oxygen vacancies formed in SrTiOs 
during the growth process followed by the spatial sep- 
aration between the charge compensating electrons and 
the vacancies^. Yoshimatsu et. al. suggested the con- 
ventional band bending mechanism based on the obser- 
vation of the core level shifts seen in photoelectron spec- 
troscopy for LaO-Ti02 interface in a soft x-ray regime^. 
The mechanism of particular interest is that of electronic 
reconstruction^ resulting in self doping of the interface 
layers with carrier density and mobility that are much 
higher than currently achieved in conventional semicon- 
ductors. The feasibility of electronic reconstruction sce- 
nario is currently under extensive investigation from ex- 
perimental as well as theoretical standpoint^ Here we 
have examined the role played by the band structure in 
electronic reconstruction of polar surfaces and interfaces. 
To avoid unnecessary complexity, we focus on simple ma- 
terials keeping all other degrees of freedom frozen. We 
chose magnesium oxide because of simple crystal struc- 
ture and high ionicity even though we do not expect MgO 
to electronically reconstruct because of the large cost in 
energy due to the large band gap. 

According to Tasker's classification scheme, the (001) 
surface of LaAlOa and (111) surface of MgO belong to 
type 3 polar termination of an ionic crystal and thus 
should not exist without a major reconstruction^. The 
instability occurs due to accumulation of the electrostatic 
potential in the film in the direction parallel to the sur- 
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FIG. 1: An artist's concept of charge neutral stoichiometric 
(a) and non-stoichiometric (b) slabs. S is a planar charge 
density (Q/A). The arrows depict electron transfer due to 
electronic reconstruction (a) and charge density redistribu- 
tion at the interface (6) due to the band structure effects. 



face normal. Each unit cell contributes the same poten- 
tial drop resulting in the electrostatic potential diverging 
with the film thickness. By virtue of its construction, 
only stoichiometric films with unalike polar terminations 
(Figl(a)) exhibit a divergent potential inside the mate- 
rial which drives the surface/interface reconstruction. In 
the non-stoichiometric structures with both film surfaces 
alike, the surface electronic structure is governed by the 
requirement of overall charge neutrality. For example, 
a defect-free (111) MgO film terminated on both sides 
with Mg is expected to be metallic in band theory be- 
cause there is not enough Oxygen in the system to fully 
oxidize the Magnesium. Given the geometry of the struc- 
ture and the Coulomb repulsion between alike charges, it 
is clear that the Mg on this surface should be 1+ just like 
the Mg on the surface of a stoichiometric (111) MgO slab 
in the case of electronic reconstruction. This approach 
kind of bypasses rather than treats the polar catastrophe 
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FIG. 2: Schematic representation of the octopolar reconstruc- 
tion of polar (111) surface of an ionic material with rock-salt 
structure. 



issue but is a convenient approximation to get insight 
into the properties of excess charge density at each inter- 
face separately that might not be feasible otherwise^' 
Note that the term non-stoichiometric is used by us in a 
rather specific sense and refers to the systems in which 
the number of positive and negative charged atomic lay- 
ers differs by one. It does not include formation of point 
defects such as vacancies nor adsorption of the foreign 
species on the surface. 

Unlike type 3, the non-polar crystal truncations are a 
single boundary problem. Here, the properties are de- 
termined by the discontinuity in the crystal potential at 
a single interface with another material or vacuum. For 
example, the Shockley surface state on (111) surfaces of 
simple metals such as Cu, Ag or Au is a result of such 
discontinuity^^. In a case of charge neutral (type 1) sur- 
faces of the ionic materials, it is the disruption of the 
Madelung potential and the change in the local symmetry 
and coordination numbers that governs a change in the 
electronic structure of a material at the surface. Some 
layered structures such as realized in TiS2 where Ti^+ 
layer is sandwiched between two S^~ layers do exhibit a 
charged surface (e.g. (001) surface) but just like type 1 
surfaces bear no dipole moment in the unit eel in a di- 
rection perpendicular to the surface. Hence, the surface 
energy is always finite. 

What are the solutions to the polar catastrophe gen- 
erated by a polar crystal truncation? As suggested by 
recent density functional calculations, a simple solution 
might prevail when the films are stoichiometric and only 
a few layers thick^^'^^. The prescription is to reduce the 
dipole moment in each unit volume by essentially con- 
verting the two adjacent planes (e.g. Mg (2+) and O 
(2-) (111) planes in MgO) into one almost charge neutral 
plane with vanishingly small dipole. Similar behavior was 
found in the calculations of thin LaAlOs overlayers on a 
SrTiOs substrate^''. These findings are consistent with 
the observation of a critical thickness in LaAlOs/SrTiOa 
heterostructures for metallic behaviour of the interface^^. 
There are, however, other plausible pathways whose con- 
tributions to the total energy of the system are indepen- 
dent of the film thickness. 

Basic electrostatic arguments suggest that the diver- 



gence of the electrostatic potential can be cured if the 
charge on the top and bottom surfaces is reduced to half 
of that of the corresponding layers in the bulk^^. Ad- 
sorption of a monolayer of foreign ions is, perhaps, the 
easiest way to achieve this without destroying the charge 
neutrality or introducing a significant structural modifi- 
cation on the surface. A notable example is a growth of 
a-MnS single crystals using iodine as a transport agent in 
chemical transport method^^. The crystals are found to 
exhibit Mn terminated (111) faces with a monolayer of I~ 
on top. Note, the crystal structure of a-MnS is rock salt 
and hence this termination should be otherwise-unstable. 
As another example, the hydroxy 1 groups have been ob- 
served on the polar surfaces of divalent oxide films^^^^. 

The surface charge can also be altered by vacancies. 
For example, the recently discovered high temperature 
superconductor BaFe2As2 has a nominal layered struc- 
ture stabilized by the extremely large electronic polar- 
izability of pnictogens^^^^. It consists of negatively 
charged (001) FeAs sandwiches alternated with the Ba 
1+ layers that generates a polar surface when the mate- 
rial is cleaved parallel to the layers. Scanning tunnelling 
microscopy studies suggest that the surface generated by 
cleaving is the Ba-terminated surface with every other Ba 
missing^^'^^. A more involved defect structure is often 
discussed in relation to (111) surfaces of the ionic ma- 
terials with a rock salt crystal structure. For example, 
the so-call octopolar reconstruction is one in which 3/4 of 
top layer atoms are removed together with 1 /4 of those of 
the layer underneath, as shown in Fig. 2, creating again a 
50% charged surface^'^. Interestingly, recent DFT stud- 
ies of octopolar reconstruction on (111) surfaces of MgO 
and NiO suggest very strong deviation from a proposed 
structure toward rather flat surfaces consisting of a sin- 
gle buckled layer of metal and oxygen ions with a formal 
charge of ±1^^. 

Electronic reconstruction is yet another compensation 
mechanism. This mechanism reminds us of the Zener 
breakdown in which electronic tunnelling between two 
interfaces occurs when the external electric field ex- 
ceeds the critical value imposed by the band gap of the 
dielectric^^. There is, however, an important difference 
between Zener theory and polar catastrophe. In the lat- 
ter case, the electric field is an intrinsic feature of an 
uncompensated film and should vanish if surface charges 
are exactly balanced with respect to charges of layers 
underneath. 

Let us consider a simple capacitor like model of a stack 
of 2N alternating charged planes with the charge density 
S and separation d as shown in Fig. 1(a). The potential 
difference between two interfaces is given by 

Vslah = ^^^{2N - l)d + (1) 

where a is the electron transfer between two interfaces 
(0 < <j < E). Here, we refer to the surface as the in- 
terface. The second term is the bulk electrostatic po- 
tential, Vbuik- When a = 0, Vsiab is simply equal to 
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FIG. 3: GGA layer projected Mg (a) and O (b) density of states for 18 bilayer thick MgO (111) slab. The zero of energy is 
at the Fermi energy, Ef- The arrows depict the lower (electrons) and upper (holes) limits in the integration about Ef used 
to calculate charge densities. Reconstructed charge distribution (c) and planar averaged total potential (d) in the slab as a 
function of distance z in direction parallel to surface normal. The residual electric field as a function of slab thickness (insert 
in (d)). A comparison between bulk and slab total density of states (e). 



twice the number of bilayers times the potential drop 
in the bulk material (e.g. potential drop between adja- 
cent Mg and O (111) planes in the unit cell of MgO). 
To put it in prospective, the 4.3 nm thick (111) film 
of MgO (18 bilayers) would experience an unsustainable 
1030.7V potential across it. It is only about a factor of 
2.4 smaller in LaAlOa (001) slab consisting of the same 
number of atomic layers (429.5V). Nevertheless, when the 
charge density on the both interfaces is reduced to a half 
(cr = ^), the potential inside the slab equals the bulk 
potential and the electric field averages out to zero sug- 
gesting that the properties of fully charge compensated 
film are, interface layers aside, representative of that in 
the bulk. While capturing the basic electrostatics, such a 
simple model is lacking some key ingredients. Neglecting 
the polarizability effects, it is clear that redistribution of 
the electronic charge between two unlike interfaces re- 
quires an excitation across the band gap. 

This brings us to the main subject of the paper: the 
role played by the band gap in electronic reconstruction 



of polar surfaces and interfaces. Let us define a quantity 
SV = Vsiab — Vbuik as a measure of the divergence in the 
electrostatic potential of the polar slab. As long as SV 
is larger than the band gap, the electrons would have to 
carry the compensating charge from one interface to an- 
other until 6V = A condition is fulfilled. Neglecting the 
band structure effects, the electron transfer is given by 
q = \Q— ^2N-i)d ^ where Q is a charge per unit area A of 
the atomic plane in the bulk. Note that A decreases the 
amount of charge transfer and only in the limit of very 
thick films {N oo) the films are fully electronically re- 
constructed {q = ^Q). Using 7.8eV band gap and assum- 
ing fully ionic charges, we estimate the amount of charge 
left on the reconstructed surface of 18 bilayers thick MgO 
(111) slab to be ±1.008e. This corresponds to a trans- 
fer of 0.992 electron between two surfaces of the slab. 
We get ±1.004e surface charge when the gap is reduced 
to DFT value of 4.5eV. The incomplete reconstruction 
leaves behind the residual electric field of A/{2N — l)d. 
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The field is greater for tiie larger band-gap materials and 
can alter the crystal structure of very thin films as was 
demonstrated recently by the DFT calculations of lat- 
tice relaxation in various polar films^^'^^'^^'^^. Note that 
there is no internal field in the non-stoichiometric case. 

Turning the discussion to the results of our DFT stud- 
ies, we note that MgO has the rock salt crystal structure 
with a lattice constant of 4.213A . It is a perfect dielec- 
tric with the band gap of 7.8eV^^'^^. Using experimental 
lattice constant, we get the band gap of 4.5 eV in agree- 
ment with previous density functional calculations^^. 
The MgO electronic structure is calculated using the 
pseudopotential DFT code SIESTA^^^ with the norm- 
conserving Troullier-Martins pseudopotentials"^^ and the 
non-polar basis set of double- C quality for s-p and single- 
( for (i-orbitals. The exchange and correlation effects 
are treated within the generalized gradient approxima- 
tion (GGA), after Perdew et al^'^ . In the slab calcula- 
tions, we adopt an ideal bulk crystal structure and pe- 
riodic slab geometry with lOOA separation between the 
slabs. The effects of periodic boundary conditions will 
be discussed latter. The stoichiometric slabs consisting 
of 6-18 MgO bilayers are used throughout this study to 
address the polar catastrophe. No lattice relaxation has 
been performed in oder to isolate the effect due to elec- 
tronic reconstruction. 

Figures 3 (a) and (b) show the layer projected density 
of states for 18 bilayers thick MgO film terminated with 
Mg layer on one side and Oxygen on the other. The pres- 
ence of the electric field is clearly evident from the shifts 
in the layer projected density of states and the slope of 
planar averaged total potential shown in Fig. 3(d). To 
find the exact amount of charge transfered between two 
surfaces, we calculate the two dimensional electron den- 
sity of the slab as a function of the distance in a direc- 
tion perpendicular to the slab. The electron and hole 
densities are calculated separately as the integrals over 
the ab-plane and over the energy from -3eV to Ep (elec- 
tron) and from Ep to +2eV (hole) as depicted by the 
arrows in Fig. 3 (a) and (b), respectively. To select the 
surface related contributions, the densities are truncated 
with a cutoff value of lO^'^e/A. The apparent difference 
in spatial distribution of hole and electron densities as 
seen in Fig. 3(c) is due to the spatial extend and orien- 
tation of the atomic orbitals these particles are in. The 
electrons occupy rather extended and isotropic Mg 35 or- 
bitals whereas the holes are in in-plane O 2p states. This 
is a result of the Hartree density-density repulsion be- 
tween the O 2p electrons which is at its minimum when 
the holes are confined to the surface. Further integra- 
tion over z reveals 0.9860±2 x 10~^ electron transfer (g). 
It is remarkable how close it is to ±1 reconstruction ex- 
pected for an ideal slab in which Mg and O are fully ionic 
with charges ±2. Reduction in slab thickness results in 
a persistent decrease in q as shown in fig. 4. A deviation 
from the estimates given above can be attributed to ac- 
tual electron density distribution and the band structure 
effects being not accounted for in the ionic model. 
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FIG. 4: Electron transfer (q) between two surfaces of the 
free standing MgO (111) film as a function of slab thickness 
obtained from DFT electron densities as described in the text. 
The solid line is a fit to the data using eq.2. 

As one can see in Fig. 3(c), the electron density on the 
Mg surface is rather extended in a direction away from 
the slab changing the electrostatic potential in the slab. 
Therefore, the expression for electron transfer needs to 
be modified. It is easy to show that, for this particular 
case, q is given by 

_ Q/2 - AeoA/D 

where D = {2N — l)d and do is the distance from Mg 
surface to the mean of the electron density in vacuum. 
Note that the equation reduces to that given above in 
a limit of do = 0. Using the 3.6V residual potential 
found in DFT calculation of 18 bilayer thick slab, we get 
Q = 1.997 ± 0.003e and do = 0.39 ± 0.03A by fitting 
eq.2 to the DFT data for electron transfer obtained for 
various slab thicknesses. 

Due to the residual electrostatic potential, the total 
density of states is broadened to, approximately, the total 
width of that in the bulk plus the gap energy. The effect 
is most evident in a case of O 2s states whose band widths 
are increased from 2eV in the bulk to about 5eV as shown 
in Fig. 3(e) suggesting that the spectroscopy of polar 
films could be very different from the bulk or films with 
non polar terminations. 

What about the periodic boundary condition? Since 
the slab is charge neutral there should be no electric field 
outside of the slab. The periodic boundary condition 
creates an artificial potential drop with the slope propor- 
tional to l/Dy where Dy is the distance between the 
slabs. This is analogous to non polar slabs with unlike 
terminations such as SrO and Ti02 surfaces in SrTiOa 
(001) slab. This is often avoided by treating the surfaces 
separately e.g. SrTiOa (001) slab terminated with SrO 
layer on both sides. It is not an option in a case of polar 
slabs. One can wonder if an artifact due to PBC has any 
influence on the electronic structure of polar films. A 
simple way to resolve PBC problem is to create a double 
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FIG. 5: MgO (111) planar averaged potential as a function 
of distance parallel to slab normal, V{z), calculated in double 
slab geometry consisting of two slabs per supercell facing each 
other. The insert shows the difference with conventional slab 
model. 



slab supercell with two slabs facing each other. In other 
words, the top surface of one slab is the same as the bot- 
tom surface of the slab above. The potential in vacuum is 
then constant which makes such a construction periodic. 
Fig. 5 shows a planar averaged potential as a function of 
distance in direction parallel to the surface normal calcu- 
lated using double slab supercell. Apart from the vacuum 
region, there is no distinction between the potentials cal- 
culated using single and double slab super cells. In both 
cases, the residual potential inside the slab is about 3.6V. 
We conclude that for sufficiently large Dy the electronic 
structure of the film calculated in a single slab geome- 
try is not affected significantly by the periodic boundary 
condition. 

In free-standing films, the main cause for incomplete 
electronic reconstruction is the films band gap. It 
is, however, a conceptual system since films are never 
grown without the substrate. Furthermore, the inter- 
faces are essential part of any heterostructure such as 
LaAlOa/SrTiOa. One can wonder if an interface with 
non-polar material can alter the electronic reconstruc- 
tion of the polar films. Neglecting the band structure 
effects, this question leads to reexamination of the en- 
ergy scales involved in the process. It is clear that, as 
long as the electronic structure on the non-polar side of 
the interface is such that it can provide states in the band 
gap of the polar one, it would have to be a part of the 
electronic reconstruction. Nevertheless, the charge redis- 
tribution within a single interface can not cure the polar 
catastrophe. Electrostatics requires that both interfaces 
participate in the same manner as in a free standing film. 

Let us now examine the electronic structure of a sto- 
ichiometric MgO (111) film on the surface of a metal, 
namely, aluminum. To keep the discussion coherent, we 
assume an ideal crystal structure and the lattice constant 
of MgO (Al is isotropically expanded). The lattice mis- 
match between Al and MgO is 3.9 percent. The structure 
is modelled by 18 bilayers of MgO and the same number 
of monolayers of Al in a single slab geometry with 100 A of 
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FIG. 6: GGA layer projected density of states for MgO (111) 
film on top of 18 monolayers of aluminum with O (left column) 
and Mg (right) contacts. Top, middle and bottom panels 
are Al surface, interface and MgO surface DOS, respectively. 
Color coding: Al is black, Mg and O are red and blue. The 
zero of energy is at the Fermi energy, Ep. 



vacuum. The interface distances are chosen according to 
the minimum in the total energy as a function of MgO- Al 
distance. For the 0-Al interface, we find that the separa- 
tion is 1.44 A whereas it is substantially larger in a case of 
Mg-Al contact, 2. 4 A. Quick inspection of layer projected 
DOS shown in Fig. 6 reveals that (i) the MgO surface 
projected DOS is very close to that in a free standing 
film (fig. 3); (ii) the residual potential is 4.0V in MgO/Al 
slab containing Mg-Al interface and it is reduced to zero 
in the case of an 0-Al contact layer; (iii) Al surface pro- 
jected density of states is not affected by the interface 
with MgO. 

It is impossible to accurately predict the band align- 
ment in the case of polar interface because the films are 
always reconstructed electronically in the calculations. 
However, the difference in 8V can be understood if the 
Fermi energy of Al was closer to the conduction band 
rather then the top of valence band of MgO. If that is 
the case, the redistribution of the charge between Mg 
surface and 0-Al interface would necessitate the energy 
scale that is less than a half of the band gap which trans- 
lates into smaller W and increases the electron transfer. 
This is exactly what we find. It is clear that having Al 
electronic states in the band gap of MgO helps to reduce 
the energy needed to reconstruct electronically because 
they provide the reservoir of the compensating charge at 
a fraction of the cost. The end result depends, however, 
on the film orientation (O vs Mg contact) and layer spac- 
ing at the interface. For example, we find that the total 



6 




1.8 2.1 2.4 

Interface Spacing (Ang.) 



FIG. 7: The residual potential in MgO (111) film on the Al 
substrate as a function of interfacial spacing (dint)- The solid 
symbols depict dint corresponding to a minimum in the total 
energy as a function of dint- 



energy of the film with Mg-Al contact is about 3.1eV 
higher as compared to 0-Al contact which is simply an- 
other evidence of much larger energy scale involved in 
the electronic reconstruction of the film. Again, we stress 
that due to electrostatic nature of the problem both sur- 
face and interface are strongly coupled via the electro- 
static potential. This drives the electron transfer from 
Oxygen surface into Mg-Al interface that, in this partic- 
ular configuration, costs almost as much as the band gap 
energy and leads to the residual potential which is nearly 
independent on the interface spacing. Fig. 7. The switch 
to another geometry allows to utilize much smaller en- 
ergy scale and reduce the residual potential to zero at 
optimal 0-Al distance. Of course, when this distance 
is large, the system is no longer coupled and the elec- 
tronic reconstruction pathway is independent on the Al 
substrate. So, the reconstruction energy scale is, again, 
the band gap of MgO. 

To further demonstrate the interplay between the 
structure and the energy in electronic reconstruc- 
tion pathways, we calculate the electronic structure 
of MgO (111) film sandwiched between two Al slabs 
(Ali2 / (MgO) 12 / Ali2 ) . The interface distances are chosen 
to be equal to that found in MgO/Al slab calculations. 
Figure 8 shows the layer projected density of states and 
the planar integrated valence electron density difference 
defined as Ap{z) = pif^M^o/Aii^) ' Pifi^) " Pm,o(^)- 
Comparing to Fig. 6, it is clear that the MgO projected 
DOS is very similar to that in MgO/O-Al/Al slab with 
the exception of the Mg surface projected DOS. In partic- 
ular, there is no shift in MgO DOS that suggests that the 
structure is fully compensated electronically. The elec- 
tron density difference plotted in Fig. 8(b) as a function of 
distance z in a direction parallel to interface normal illus- 
trates the role of Al band structure in the reconstruction 
pathway. As in other structures considered, the electron 
transfer is across the MgO slab but it is, to a large ex- 
tent, Al density that participate in the process. This 
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FIG. 8: Left panel - GGA layer projected DOS for 
Ali2(/MgO)i2/Ali2 slab. The zero of energy is at the Fermi 
energy. Middle panel - valence electron difference Ap as a 
function of distance z in direction parallel to interface normal 
(Apiz) = ptiM90iAi{^) - Pm\^) - PMgoiz)). The oscil- 
lations inside MgO fragment of the structure are due to in- 
complete electronic reconstruction in free standing MgO (111) 
film. The oscillations at the interfaces are due to the change 
in reconstructed electron (hole) density distribution as com- 
pared to single MgO slab. Right panel - an illustration of the 
change in spatial distribution of compensating electron and 
hole densities due to interface with Al. The large oscillations 
at the interfaces in (b) are primarily due to readjustment of 
the compensating densities as compared free standing film. 



reduces the double interface energy {Edint) to 1.92eV as 
compared to the slab surface energy (Esse) of 6.02eV in 
a free standing MgO. The Edint is given by 



E^ 



_ r^Al/MgO/Al 
dint — -t^tot 



-NE, 



Al bulk 



where N (M) is the number of formula units in Al (MgO) 
slab; Etot is the respective total energy; i^ surface 

energy calculated for 12 monolayer thick Al (111) slab us- 
ing MgO lattice constant (0.42eV). Correspondingly, the 

1 -1 r -771 T-iMqO slab t\ /r -i-iMaO bulk 

slab surface energy is Esse = E^^^ — ME^^^ 

In conclusion, we demonstrate that the shifts in the 
DOS of polar systems are extremely sensitive to tiny de- 
viations in electron transfer involved in the electronic re- 
construction mechanism. They are not the signatures of 
electronic reconstruction and only tell us about its in- 
completeness. It is a result of an effective energy scale 
which is required in order for charge to be transferred be- 
tween the oppositely charged interfaces. In a free stand- 
ing film, this energy scale is set by the band gap modified 
to some degree to incorporate the band structure effects. 
At interfaces, it is determined by the band alignment 
and it is always the smallest energy possible helping the 
electronic reconstruction mechanism to be realized more 
readily in heterostructures as compared to free standing 
films. Due to high sensitivity of the inner potential vari- 
ation to the interface charges, one can expect that the 
interface imperfections such as vacancies or adsorption 
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of foreign species can reduce significantly the residual 
electrostatic potential even though the electron transfer 
energy scale of pure material is still rather large^^. This 
can, in principle, reconcile theory and the results of x- 
ray photoelectron spectroscopy measurements in which 
no inner potential variation was found^^'^^'^^'^"^ Fur- 
thermore, bulk impurity states energetically distributed 



across the band gap can somewhat modify the charge 
distribution at the interface in a fashion similar to con- 
ventional semiconductor P-N junctions. 
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